Transverse-Momentum Resummation for Slepton-Pair Production at the LHC 
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We perform a first precision calculation of the transverse-momentum (qr) distribution of slepton 
pair and slepton-sneutrino associated production at the CERN Large Hadron Collider (LHC). We 
implement soft-gluon resummation at the next-to-leading logarithmic (NLL) level and consistently 
match the obtained result to the pure fixed-order perturbative result at leading order (LO) in 
the QCD coupling constant, i.e. 0{as). We give numerical predictions for fifi and fiiy* + t^Ut 
production, also implementing recent parameterizations of non-perturbative effects. The results 
show a relevant contribution of resummation both in the small and intermediate gr-regions and 
little dependence on unphysical scales and non-perturbative contributions. 



INTRODUCTION 

The Minimal Supersymmetric Standard Model 
(MSSM) [1, 2] is one of the most promising extensions 
of the Standard Model (SM) of particle physics. It 
postulates a symmetry between fermionic and bosonic 
degrees of freedom in nature and predicts the existence 
of a fermionic (bosonic) supersymmetric (SUSY) partner 
for each bosonic (fermionic) SM particle. It provides a 
qualitative understanding of various phenomena in par- 
ticle physics, as it stabilizes the gap between the Planck 
and the electroweak scale [3], leads to gauge coupling 
unification in a straightforward way [4], and includes 
the lightest supersymmetric particle as a dark matter 
candidate [5]. Therefore the search for supersymmetric 
particles is one of the main topics in the experimental 
program of present (Fermilab Tevatron) and future 
(CERN LHC) hadron coUiders. 

SUSY must be broken at low energy, since spin part- 
ners of the SM particles have not yet been observed. As 
a consequence, the squarks, sleptons, charginos, neutrali- 
nos and gluino of the MSSM must be massive in compar- 
ison to their SM counterparts. The LHC will perform a 
conclusive search covering a wide range of masses up to 
the TeV scale. Total production cross sections for SUSY 
particles at hadron colliders have been extensively stud- 
ied in the past at leading order (LO) [6-8] and also at 
next-to-leading order (NLO) of perturbative QCD [9-15]. 

We focus our attention on slepton pair (slepton- 
sneutrino associated) production at the LHC through the 
neutral (charged) current Drell-Yan (DY) type processes 
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Due to their purely electroweak couplings, sleptons are 
among the lightest SUSY particles in many SUSY- 
breaking scenarios [16] . Sleptons and sneutrinos often de- 
cay directly into the stable lightest SUSY particle (light- 
est neutralino in mSUGRA models or gravitino in GMSB 
models) plus the corresponding SM partner (lepton or 



neutrino). As a result, the slepton signal at hadron col- 
liders will consist in a highly energetic lepton pair, which 
will be easily detectable, and associated missing energy. 

In this Letter, we study the transverse- momentum {qx) 
distribution of the slepton pair. Since in hadronic colli- 
sions the longitudinal momentum balance is unknown, a 
precise knowledge of the (/^-balance is of vital importance 
for the discovery of SUSY particles. In the case of slep- 
tons, the Cambridge (s)transverse mass mT2 proves to be 
particularly useful for the reconstruction of their masses 
[17] and determination of their spin [18], the two key fea- 
tures that distinguish them from SM leptons produced 
mainly in WW or tt decays [19, 20]. Furthermore, both 
detector kinematical acceptance and efficiency depend, 
of course, on qT. 

When studying the gr-distribution of a slepton pair 
produced with invariant mass M in a hadronic collision, 
it is appropriate to separate the large-gr and small-gr 
regions. In the large-gr region {qr > M) the use of 
fixed-order perturbation theory is fully justified, since the 
perturbative series is controlled by a small expansion pa- 
rameter, a,{M^). The QCD [14] and full SUSY-QCD 
[15] corrections for slepton pair production are known 
to increase the hadronic cross sections by about 25% at 
the Tevatron and 35% at the LHC, thus extending their 
discovery reaches by several tens of GeV. Recently the 
LO calculation for slepton pair production has been ex- 
tended to include mixing between left- and right-handed 
sfermions and longitudinal beam polarization [21]. 

The bulk of the events will be produced in the small-gr 
region, where the coefficients of the perturbative expan- 
sion in as(M^) are enhanced by powers of large logarith- 
mic terms, ln(M^/g|.). As a consequence, results based 
on fixed-order calculations diverge as — > 0, and the 
convergence of the perturbative series is spoiled. These 
logarithms are due to multiple soft-gluon emission from 
the initial state and have to be systematically resummed 
to all orders in a., in order to obtain reliable perturbative 
predictions. The method to perform all-order soft-gluon 
resummation at small qx is well known [22-32]. The re- 
summation of leading logs was first performed in [22]. 
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It was shown in [23] that the rcsumniation procedure 
is most naturally performed using the impact-parameter 
(6) formalism, where b is the variable conjugate to qt 
through a Fourier transformation, to allow the kinemat- 
ics of multiple gluon emission to factorize. In the special 
case of DY lepton pair or electroweak boson production, 
&-space resummation was performed at ncxt-to-lcading 
level in [25], an all-order resummation formalism was 
developed in [29], and the next- to-next-to- leading order 
terms have been calculated in [30]. 

At intermediate qx the resummed result has to be con- 
sistently matched with fixed-order perturbation theory in 
order to obtain predictions with uniform theoretical ac- 
curacy over the entire range of transverse momenta. 

In this work we implement the formalism proposed in 
[31, 32] and compute the (^T-distribution of a slcpton pair 
produced at the LHC by combining NLL resummation at 
small qt and LO {0{as)) perturbation theory at large qt- 



Qt-RESUMMATION at the NLL LEVEL 

The partonic cross section for DY slepton pair produc- 
tion can be written as 
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where a, b label the partons wliic:li take part in the hard 
process. The resummed contribution can be written as 



(3) 



^^(ros.) 

-^jj^0^iqT,M,s-,as{fJ,%),f4i,fj,%) = 
f°° b 

— db- Joibqr) Wabib, M, s; as{i4i), lA^ Mf) , 

where Jo (a;) is the 0*''-order Bessel function, {hf) 
is the renormalization (factorization) scale, and s is the 
partonic center-of-mass (CM) energy. 

The perturbative function W embodies the all-order 
dependence on the large logarithms \n.{Al'^b'^). They cor- 
respond, in the conjugate space, to the previously men- 
tioned terms, \n{M'^ /q^), that spoil the convergence of 
the perturbative series at small qT (large b). Perform- 
ing a Mellin transformation with respect to the variable 
z = M"^ /s at fixed M, we can define the AT-moments Wn 
of W and express them in an exponential form 



Wjv(6,M;a4/4),/i|„/4) = 

Hn {M,as{iil);Mytxl,Myfil,MyQ^) x 

eMGN{aM),L;Myfil,MyQ')} , 
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where constant (i.e. finite as qx 0) and logarithmically 
divergent terms are factorized into the functions Hn and 
Qn, respectively. This factorization implies some degree 
of arbitrariness, and the scale Q is introduced to param- 
eterize this uncertainty. As in the case of and hf, one 



should sot Q = M and estimate the uncertainty from un- 
calculated subleading logarithmic corrections by varying 
Q around this central value. 

The function Tijy docs not depend on the impact pa- 
rameter b and, therefore, it contains all the perturbative 
terms that behave as constants in the limit 6 — > oo. In 
addition it contains the whole process dependence as well 
as factorization scale and scheme dependence. Its expan- 
sion in powers of as gives 



(5) 



aW(a„M) 1 + ^7^)^^(MVM'fl,MV/x2^,M7Q2) 



/(I) 



h<;'(m7pJ„mVm5-,m79') + .^^ 



where ct^°^ is the lowest-order partonic cross section for 
the hard-scattering process. The coefficient Ji^^ splits 
into a process-independent flavor off-diagonal contribu- 
tion and a process-dependent [33, 34] flavor diagonal con- 
tribution. The general expression for , needed to per- 
form a NLL analysis, can be found in [32]. The second 
order coefficient n^jq has not yet been computed. 

The exponent Qj^ includes all the terms that are log- 
arithmically divergent when 6 ^ oo (gr 0) (i.e. pro- 
portional to Z/ = ln(Q262/6o), where 6o = 2e~'>'® and 7^ 
is the Eulcr number) . However these terms become large 
both for small and large 6-values, introducing unjustified 
large contributions also at large qx- It is useful to intro- 
duce a modified expression of the expansion parameter. 



L = In 



2l2 



+ 1 



(6) 



The variables L and L are equivalent at small q^ , but dif- 
fer at intermediate and large qx, avoiding the unwanted 
resummation contributions since Z — > for ^ 1 and 
allowing us to recover the corresponding fixed-order total 
cross section after integration over qx- 

The exponent Q can be systematically expanded as 

ejv(a«,L;MVM|s,MVQ2)= (7) 
^5i'^(a,L;MV/.LMVQ') + ..., 

TT 

where the term Lg^^^ collects the leading logarithmic 
(LL) contributions, the function ^^^^ resums the next-to- 
leading logarithmic (NLL) contributions, a^g'-'^-' controls 
the next-to-next-to-leading logarithmic (NNLL) terms, 
and so forth. The explicit expressions for the gi functions 
arc given in [32] in terms of the universal perturbative co- 
efficients 1 ^q^N^ ^'<h^N ' Sincc wc want to 
perform the resummation at the next-to-leading level, we 
need to implement the LL fimction Lg^^\ which depends 
on the coefficient A^q \ and the NLL function Lg'^'^^ , which 
depends on A^q^ and -B^^^. 
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The second term [daf^ '^ /dM'^dql) in Eq. (2) is free 
of divergent contributions and can be computed by fixed- 
order truncation of the perturbative series. In order to be 
consistently matched with the resummed contribution at 
intermediate qr {ffr — M) , this term should be evaluated 
starting from the usual perturbative truncation of the 
partonic cross section and subtracting from it the expan- 
sion of the resummed part at the same perturbative or- 
der. This matching procedure between small- and large- 
qt regions prevents double-counting (or neglecting) of 
perturbative contributions and guarantees a uniform the- 
oretical accuracy over the entire transverse-momentum 
spectrum. Since the fixed-order cross section for slep- 
ton production at non-vanishing transverse-momentum 
is known at LO (slepton pair -I- one parton) [14, 15], we 
can only consistently perform a NLL-I-LO matching. 

The above formalism refers to a purely perturba- 
tive framework. Nonetheless it is known [27] that the 
transverse-momentum distribution is affected by non- 
perturbative (NP) effects which become important in 
the large-& region. In the case of electroweak boson 
production, these contributions are usually parameter- 
ized by multiplying the function W in Eq. (3) by a NP 
form factor T^^(h) [35-39], whose coefficients are ob- 
tained through global fits to DY data. We include in 
our analysis three different parameterizations of NP ef- 
fects corresponding to three different choices of the form 
factor: the Ladinsky-Yuan (LY-G) [36], Brock-Landry- 
Nadolsky-Yuan (BLNY) [38], and the recent Konychev- 
Nadolsky (KN) [39] form factor. 



SLEPTON PAIR PRODUCTION AT THE LHC 

In this Section we present quantitative results for the 
^T-spectrum of slepton pair (slepton-sneutrino associ- 
ated) production at NLL+LO accuracy at the LHC col- 
lider. We focus our study on the lightest slepton mass 
eigenstate fi and thus we consider the processes 



p p ^^,Z^X,T.* at the LHC 



qq' ^ Tiv*,flvr- 



(8) 



We use the MRST (2004) NLO set of parton distribution 
functions [40] and as evaluated at two-loop accuracy. We 
fix the resummation scale Q equal to the invariant mass 
AI of the slepton (slepton-sneutrino) pair and we allow 
ji — jip — IJ.R to vary between M/2 and 2M to esti- 
mate the perturbative uncertainty. We also integrate Eq. 
(2) with respect to M^, taking as lower limit the energy 
threshold for Tifj*(fi£',-) production and as upper limit 
the hadronic energy (75=14 TeV at the LHC). 

In the case of fifl production (neutral current pro- 
cess, see Fig. 1), we choose the SPS7 mSUGRA bench- 
mark point [16] which gives, after the renormalization 
group (RG) evolution of the SUSY-breaking parameters 
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FIG. 1: Differential cross section for the process pp fif]*. 
NLL-I-LO matched result, LO result, asymptotic expansion of 
the resummation formula and A-parameter are shown. 



performed by the SUSPECT computer program [41], a 
light fi of mass m^j = 114 GeV. 

In the case of fii)* +flDT production (charged current 
process, see Fig. 2), we use instead the SPSl mSUGRA 
benchmark point which gives a light fi of mass ^ = 136 
GeV as well as a light Vr of mass m^^ — 196 GeV. 

In both cases we plot the LO result (dashed line), the 
expansion of the resummation formula at LO (dotted 
line), the total NLL-I-LO matched result (solid line), the 
uncertainty band from scale variation, and the quantity 



A = 



^^(rcs.+NP)(^ = M) - dcr('-'="-)(/^ = M) 

da'^^'^^'-Hp = M) ■ 



(9) 



The parameter A gives thus an estimate of the contri- 
butions from the different NP parametrizations (LY-G, 
BLNY, KN) that we included in the resummed formula. 

We can see that the LO result diverges to -|-oo, as ex- 
pected, for both processes as qx 0, and the asymp- 
totic expansion of the resummation formula at LO is 
in very good agreement with LO both at small and in- 
termediate values of qx- The effect of resummation is 
clearly visible at small and intermediate values of qx , the 
resummation-improved result being nearly 39% (36%) 
higher at qx — 50 GeV than the pure fixed order result 
in the neutral (charged) current case. When integrated 
over qr, the former leads to a total cross section of 66.8 
fb (12.9 fb) in good agreement (within 3.5%) with the 
QCD-corrected total cross section at 0[as) [15]. 

The scale dependence is clearly improved in both cases 
with respect to the pure fixed-order calculations. In the 
small and intermediate gr-region (up to 100 GeV) the 
effect of scale variation is 10% for the LO result, while it 
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pp ^ tJ V, + T, V* at the LHC 
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FIG. 2; Same as Fig. 1 for the process pp — > fiu* + t^Vt- 

is always less than 5% for the NLL+LO curve. Finally, 
non-perturbative contributions are under good control. 
Their effect is always less than 5% for qt > 5 GeV and 
thus considerably smaller than resummation effects. 

CONCLUSIONS 

In this Letter, a first precision calculation of the qt- 
spectruni for SUSY particle production at hadron collid- 
ers has been performed by applying the gy-resummation 
formalism at the NLL+LO level to slepton pair and 
slepton-sneutrino associated production at the LHC. The 
numerical results show the importance of resummed con- 
tributions at small and intermediate values of gr, both 
enhancing the pure fixed-order result and reducing the 
scale uncertainty. 
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